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Abstract 

^L ■ We study the stability of supersymmetry breaking vacuum in the presence of cosmic strings 

^vq . arising in the messenger sector. For certain ranges of the couphngs, the desired supersymmetry 
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breaking vacua become unstable against decay into phenomenologically unacceptable vacua. This 
sets constraints on the range of allowed values of the coupling constants appearing in the models 
and more generally on the chosen dynamics of gauge symmetry breaking. 
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I. INTRODUCTION 



Supersymmetry is a powerful principle which brings control in Quantum Field Theory. In 
particular its non-renormalization properties[l||2] and its potential for solving the hierarchy 
problem |3|, Ij, |5|, la] make it appeahng as an ingredient of phenomenology, enabling a descrip- 
tion in terms of weakly coupled degrees of freedom. In practice however, the implementation 
of the principle requires breaking supersymmetry in a hidden sector and communicating the 
effects to observed phenomenology indirectly; for reviews see |7|, [Sj. The issue is further com- 
plicated by the fact that breaking supersymmetry is not itself a generic possibility. It was 
elucidated by |9| that presence of R parity and its breakdown can play a crucial diagnostic 
role for understanding generic possibilities of breakdown. A persistent problem of supersym- 
metry breakdown however is that such vacua may be only local minima and have the danger 
of relaxing to a true minimum where supersymmetry is restored. This unpleasant possibil- 
ity is however easy to avoid if the tunneling rate to the true vacuum can be made much 
longer than the known age of the Universe. New insights into consistency of supersymmetry 
breakdown in metastable vacua have been obtained recently in 10|, lll[ in the context of su- 
persymmetric QCD. In such models several possibilities exist, e.g., supersymmetry breaking 
vacuum need not be unique but could be degenerate with several supersymmetry preserving 



vacua. The issues involved in phenomenologica 
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191]. Applications to cosmology appear 



The question of whether the metastable vacua are cosmologica! 
tunneling for all ranges of couplings received early attention in 



ly stable against quantum 
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in the context 



of Gauge Mediated Supersymmetry Breaking (GMSB) (see for instance the reviews 25] 



261]). In these models supersymmetry is dynamically broken in a hidden sector at a scale A^ 



and communicated to the standard model through a "messenger" sector at TeV scale. The 
messenger sector has suitable interactions with both the hidden and visible sectors. The true 
vacuum in some of these models is supersymmetric and color breaking. Supersymmetry is 
broken and colour is preserved only in metastable vacua. 

In this paper we advance the point that it is not sufficient to study the stability of 
supersymmetry breaking vacua in their translationally invariant i.e., spatially homogeneous 
avatar. Generically topological defects such as cosmic strings or monopoles may obtain in 



the course of implementation of any particular scheme. Such defects could arise in the early 
universe [23]. It is then necessary to take account of their presence when studying stability 
issues. A class of topological defects can exist which can nucleate the formation of the true 
vacuum in such a way that the exponential suppression inherent in tunneling phenomenon 
no longer obtains. Specifically, when a cosmic string is present in a false vacuum, such a 



local minimum can be rendered unstable against decay to a vacuum of lower energy 28| 29 1. 



This process was studied by one of us 30|] in the context of phase transitions in Grand 
Unified Theory (GUT) models. It was shown that generically, the presence of the cosmic 
string entailed the consequence that false vacuum would "roll-over" smoothly to the true 
vacuum without recourse to quantum tunneling. In this way, a putative first-order phase 
transition becomes second order, with important cosmological implications. 

The same process is applicable here in the context of gauge mediated supersymmetry 
breaking models. It is possible to establish the required results numerically where formal 
methods provide a suggestive answer. We study a model of supersymmetry breaking which 
contains two classes of vacua, color breaking and color preserving. The color preserving 
metastable vacuum becomes parametrically unstable in the presence of cosmic strings. This 
puts constraints on the model of messenger sector being used and may make some models 
unviable. 

This paper is organized as follows. The next section [ITl provides a brief review of the class 
of models we have studied here. Section IIIII discusses the vortex (strings) ansatz possible in 
the relevant vacua and the equations of motion to be solved. The results of the numerical 
solutions are described in section llVi A semi-analytic approach explaining the numerical 
results is discussed in section |Vl followed by concluding remarks in section IVI[ 

II. A MODEL FOR MESSENGER SECTOR 



In this section, we summarize the models of 3l|, |32| as discussed in 22] . Supersymmetry 
is broken dynamically in the hidden sector. The messenger sector shares a symmetry Gm 
with the hidden sector. It is sufficient for Gm to be global, but we consider the possibility 
of Gm being local, permitting cosmologically viable strings. The case of global Gm is also 
commented upon. The messenger sector is responsible for transmitting supersymmetry 
breaking from the hidden sector to the visible sector. This messenger sector consists of the 



following: a gauge singlet, S, a pair of messenger quarks, q and q, belonging to the 3 and 
3 representation of the SU{3)c color group, and a pair of chiral superfields, A^ and P in 
vector-like representations of Gm- The possibility of anomalies in the hidden sector forces 
the introduction of extra superfields Ei to cancel the anomalies. However, in this paper, we 
consider only 'Minimal' models which do not include these extra fields. The most general 



superpotential for the messenger sector is given by 22 1 

Wmes = f^Sqq + ^S^ + XiPNS + Wi{P, N, E,). (1) 

The three coupling constants k, A, and Ai are all positive by a suitable phase definition of 
the fields. In the case where Gm is U{1), if e is the U{1) gauge coupling and P, N, and E^ 
have charges +1, —1, and yi, the effective potential for the charged scalars after integrating 



out the hidden sector is given by [22 1 



VsB = M^{\P\^ - |iV|2 + y,|E,n + M^dPp + |iV|2 + y2|^^|2) + , _ _ (2) 

in which the relation between the mass parameters Mi, i = 1,2 and the supersymmetry 
breaking scale A^ is 



Mf = ^^Kijj^r- (3) 



;2 

The factors q {i = 1,2) are of order unity, with their signs dependent on the content of 
the dynamical symmetry breaking sector. There are higher-dimensional terms coming from 
more loops and these are what the ellipsis in the expression for Vsb stand for. In addition, 
there are contributions to the scalar potential coming from the f/(l) D-term and F-terms. 
In what follows, we will consider the set of minimal models which satisfy ^^ = ^^ = 0. 
This set includes models in which there are no E fields, models in which Wi = 0, and models 
in which the E fields do not couple to P and N. After including the U{1) D-terms and F- 
terms and setting Wi = 0, we consider an explicit model in which (Mi = 0), and the scalar 
potential of the messenger sector becomes [22|] 

Vmes = j{\P\'-\N\'y + iM'2+^l\S\')i\P\' + \N\') 

+ K''\S\\\q\^ + \q\') + \Kqq + \S^ + \,PN\\ (4) 

We shall be interested in the case when M| < which signals the breakdown of the U{1). 
This makes Vmes unbounded from below, but in the total potential, higher-order terms in 



VsB result in a deep global ininimum far away in field space in which the visible sector is 
supersymmetric. For a viable local minimum, one must set g = g = in the expression for 
Vmes and in order to have S y^ simultaneously, we must have 

A > Ai. (5) 

The local minima lie at g = g = and 



\S\' = -M'.TTi7^^^^ (7) 



1-Ai/A 
A?(2 - Ai/A) 

Arg{PNS*^) = vr (8) 



Along with the condition A > Ai, stability of the local minima require the following relations 
between the couplings [22|: 

Xl < 2Xe^ (9) 

A, < ^ (10) 



Condition fITOl) is also mentioned in 33j wherein it is noted that there is still a possibility of 
color breaking vacua when (ITOl) holds. 

The scalar potential of the messenger sector given in equation (j3]) contains two important 
types of vacua. One of these, which will be denoted as \Vi) later in the paper, is the 
supersymmetry breaking local minimum described above. This vacuum has (S) ^ and 
(q) = (q) = 0) which means that SUSY is broken while the color gauge group is unbroken. 
The other vacuum, which shall be denoted IV2), also has (5*) 7^ but the fields g and g get 
non-zero VEVs. This means that | V2) is a phenomenologically undesirable minimum in which 
SUSY is still broken but the color gauge group is also broken. If cosmic strings are supported, 
these vacua are modified to what will be denoted as \Vi "'^^') and | Vg ""^■^) respectively. 
It will be shown in a later section that \Vi "'^^''^ becomes parametrically unstable towards 
decay into [Vg ""'^ ). The next section describes the available string solutions. 

III. THE STRING ANSATZ 

With the scalar potential for the messenger sector as given in (jlj), the full Lagrangian for 
all fields (P, A^, S, g, g) can be written down. The fact that string solutions are cylindrically 



symmetric makes it convenient to use cylindrical coordinates. We ignore the z-dependence 
and look for time-independent solutions. In a local minimum, \P\ = \N\ from ([6]) and (5*) is 
as given in ([7]) . The ansatz functions for the scalar fields are a simple generalization of the 



Abrikosov-Nielsen-Olesen string 3J] 35|, (see also 36l ]) 



P = r]f{r)e'' (11) 

N = r]f{r)e-'^ (12) 

S = r]g{r)e'''/^ (13) 

q = q = 7]h{r)e^l'^ (14) 

The value of t] is given by the right hand side of equation ([6]) Note that the winding number 
for the P and N fields has been chosen equal and opposite. Single valuedness requires that 
the phase complete a circuit of 2ti at infinity. This requires the phase to be an integer 
multiple of cylindrical angle 9. Along with the choice of phase for S as in (fT3l) . this winding 
ensures the condition (jHj). This is the expected minimal energy configuration. The phases 
for g, and q satisfy the expectation that their product remains negative. At the core of the 
vortex, the fields P and N must vanish, and hence /(O) = 0. At infinity, P and N approach 
their vacuum expectation values which are denoted by r/, and hence f{oo) = 1. Similarly, 
g{r) approaches the value given in (jTj) at infinity but can be non-zero in the core. The 
value of h{r) can be non-zero in the core also, but its value at infinity is of importance 
in determining whether the resulting vacuum is phenomenologically viable or not. This is 
because the vacuum with h and hence q and q non-zero is color-breaking. 

A vortex solution in P and N also includes a coupling of the fields to a gauge field A{r). 
The gauge field behavior is described by a function a{r) defined through 

Ae{r) = -a{r) (15) 

er 

Aq = Ar = A, = (16) 

where for continuity of Ag{r), we have a(0) = 0, and e is the unit of the abelian Gm charge. 
At infinity, Ag{r) is pure gauge and goes as 1/r and hence a(oo) = 1. The Lagrangian for 
the system can be written as: 

L = li^,Sr + ^{^,qy + ^i^,qr + ^\i^, + ^eA,)P\' + ^\{^, + ^eA,)N\'-^^^^ 

(17) 



where Vmes is as given in (jl]). After substituting equations (fTT]-[T5!) into the Lagrangian, the 
equations of motion can be written in terms of /(r), a{r), g{r), and h{r). They are a set of 
coupled second-order and non-hnear ordinary differential equations given by: 
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(18) 

(19) 
(20) 
(21) 



dr"^ r dr 

The fields P and N in which we have setup string configurations have a f/(l) symmetry. 
As is well known, existence of cosmic string solutions requires that the space of equivalent 
vacua be necessarily multiply connected. Generically in a larger gauge group such as SO{N), 
Spin{N) or SU{N), cosmic strings are in principle possible for P and A^ in any representation 
other than the vector or the fundamental. Thus, the conclusions of this paper have relevance 
to more general cases. 



IV. TESTING THE VACUA 



A. Homogeneous vacuum 



We begin with calculating the field values at infinity which can be obtained easily from 
equations (fT8] - [2T]) . At r = oo, all the derivative terms vanish due to translational invariance 
and so do the terms containing (1/r^). The equations get reduced to the following simple 
set of simultaneous equations: 



Ml 

|2 



rj^ 



-Ai(A-Ai)/-Aifi;/i2 



2e'f'{a-l) 



2 p3 



f + Kf 







[2Ai(Ai - X)f + 2k{k + X)h^] g + 2AV = 
[k{k + A)/ - nXif] h + k^h^ = 



(22) 

(23) 
(24) 
(25) 



One must note that values for three couplings (k. A, Ai) must be given for each unique set of 
vacuum solutions. The set of vacuum solutions consists of four solutions. One is the trivial 



solution in which all fields vanish. Another is one in which / and a are non-zero but SUSY 
breaking is not being communicated since g and h are both zero. 

The other two translation invariant minima shall be denoted |V^i) and IV2). Their field 
content is as follows: 

\Vi): f = a = l g^O h = (26) 

\V2):f<l a = l g^O h^Q (27) 

\Vi) is the desired supersymmetry breaking and color preserving local minimum described 
in section 2. The value of g is exactly as required by equation ([7]), and it is a function of 
the couplings. IV2) is an undesirable minimum in which h and hence q and q get VEVs. 
Furthermore, what is interesting is that for all ranges of couplings, its energy is always lower 
than that of \Vi). 

(y2\V^es\V2) < {Vi\Vmes\Vl) (28) 



The expression for Vmes in terms of {f,g,h) can be written as 



''mes V 



2f ( A?/ - ^ ) + 2K'g'h' + {X,f - A/ - ^h'^' 



(29) 



It is a simple exercise to check that for each set of solutions, Vmes is always smaller for IV2) 
compared to \Vi). 

B. Vacuum with gauge string 

A numerical strategy exists for solving equations (TTT] - [T5l) . This consists of using relax- 
ation techniques, after discretizing the equations of motion on a grid converting them to a 
set of coupled polynomial equations. In conformity with physical expectations the initial 
guess is one with / = a = in the core and with g and h arbitrary. At infinity, the trial is 
chosen to be either settling into \Vi) or IV2). The exact values of {f,a,g,h) in \Vi) and IV2) 
can be found by first solving equations (1221- [251). We denote by ^y^ """^ ) the static state 
of the system in which the cosmic string configuration approaches \Vi) (upto a ^-dependent 
phase) as \r\ — >■ 00. Similarly, [Vg ""■^^) denotes a static string solution which asymptotes 
to I V2). An example of the solutions in a case where both | V^ ^*"-9^^ g^^^j |y^w ""9j^ ^^^ stable 
is shown in figures [T]-|ll 

It may happen that for some values of the couplings, a solution with the expected field 
values at infinity cannot be obtained. When this happens, that solution with its asymptotic 
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FIG. 1: f(r) and a(r) for ^y^^ "^^'^ with boundary conditions as in eqn. (j26p . and k = 1.3, A 
1.45, Ai = 0.65. Here, foe = 1 and a^c = 1- 



vacuum values is not tenable. We have explored the stability of ^V^ """^ ) and \V2 ""'^) 
over a range of the couplings. For the entire range of couplings studied, the string solution 
m ""^ ) is always obtained. However, no stable static solution \Vi """^ ) can be obtained 
for certain ranges which are indicated in table [H Note specifically that for k = 1 and A = 1.1, 
values of Ai < 0.1 which are of phenomenological interest are immediately ruled out. It is 
for these ranges of the couplings that the corresponding metastable vacua are no longer local 
minima due to the presence of cosmic strings. 

In the context of the early universe, the parameters in the effective action are as a rule 
temperature dependent. The above results then imply that we may start with a phase 
wherein a state of the type |V^/** ""■b)'^ may be a local minimum but with reduction of tem- 
perature, the string solution may not exist in the sense argued above. How this comes about 
in the present example is taken up at the end of sec. |Vl The subsequent dynamics can be 
simulated by restoring time dependence in equations flT8] - [2T|) . It was shown in 30| that 
the string configuration rendered parametrically unstable undergoes a real time evolution 
into a possible stable string configuration with a modified vacuum at infinity. This can be 
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FIG. 2: g(r) and h(r) for ^y^^ "^^'^ with boundary conditions as in eqn. (j26p . and k = 1.3, A = 
1.45, Ai = 0.65. Here, g^ = 0.497317 and h^c = 0. 

referred to as a "roll-over" process, a prompt, semi-classical evolution rather than tunneling 
by spontaneous formation of bubbles. 



C. Global cosmic strings 

Our demonstration though specific to messenger group of this example may have relevance 
to a more general setting where some sector of the theory possesses global symmetries which 
are broken in the desirable vacuum. In this case, the strings arising are global cosmic 
strings. In our example this corresponds to setting a{r) identically to zero. We have looked 
for solutions without the a{r) field for some ranges of the couplings. The results are shown in 
table ini It is found that this condition rules out |Vj^^'^ "^9>'j ^^j- g^ larger range of couplings as 
compared to the case in which a{r) is present. The results for m ^*"'^-^) remain unchanged. 

A global cosmic string network is not desirable for stable cosmology 27|] . This is because 
the energy per unit length of a global string is divergent and can dominate the energy 
density of the Universe. However the divergence of the energy with distance from the string 
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FIG. 3: f(r) and a(r) for ^Y^'"''"^^)^ with boundary conditions as in eqn. (j27p . and k = 1.3, A 
1.45, Ai = 0.65. Here, foo = 0.999405 and Ooo = 1- 



is logarithmic and a transient network of such strings cannot be ruled out. We expect such 
a network to eventually relax to a homogeneous vacuum. Even a transient network however, 
may be capable of destabilizing a particular local minimum according to the discussion of 
this subsection. 



V. STABILITY ANALYSIS OF |y/^*""s)) 

In the previous section, the stability of \V^^ ""'^) has been studied numerically. As 
illustrated in tabled there are certain domains in the parameter space of n, Ai, and A in 
which \Vi ""^■^) is not admissible. The non-existence of such solutions was then understood 
as an unavailability of iv^^*^* ""^-'^ for those parameter values. This numerical result can be 
confirmed by a semi-analytic treatment following an approach discussed in 28 1. 

A preliminary analysis can be made for the translation invariant vacuum to study the 
effect of changing A. Referring to equation fl2^ . the effective squared mass ^g// ^^^ the h 
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FIG. 4: g(r) and h(r) for |V^^* ""^') with boundary conditions as in eqn. (j27p . and k = 1.3, A 
1.45, Ai = 0.65. Here, goo = 0.481526 and h^o = 0.094434. 



field is given by 



m. 



// 



2g^{K'^ + kX) -2KXif 



(30) 



Now if K and Ai are held fixed, when A is large enough, this quantity is positive and h = 
is admissible as a local minimum. As A reduces, there comes a point when m'^rr becomes 
negative. This corresponds to the regions of instability shown in table [H For example, when 
K = 1.3 and Ai = 0.65, reducing A below 1.43 drives its effective squared mass negative 
making h = unstable. 

We now proceed along the lines of 28| to analyze the linear stability of \Vi ""'^) by 



adding a small time-dependent term to the time-independent numerical solution. If small 
oscillations around this solution have only real frequencies, the configuration is considered 
stable. For a nonavailability of | y^"-^ ""^•' ^ ^ at least one of the modes of oscillation must 
possess an imaginary frequency. The stability analysis for \Vi ""^■') is greatly simplified by 
the following observation. The fields / and g do not differ significantly in states \Vi """^'^ 
and \V2 ""'^ ) as seen from figures [T]-|ll Only h differs significantly and the possible time- 
dependence of the background fields / and g can be ignored. The equation of motion for h 
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TABLE I: Testing the existence of local minima viable in the presence of a cosmic string. Tick 
marks indicate when solutions are found and crosses indicate when solutions cannot be obtained. 
Vacua of type 1, eqn. (|26|) . become parametrically disallowed due to the presence of a cosmic 
string. 
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TABLE II: Sample results for global strings, i.e., no a{r) field present. These are to be compared 
to rows 5 - 9 of table HI 



with time-dependence restored is 






and we now write 



h{r,t) = h{r) +p(r)e 



iujt 



(31) 



(32) 



13 



Here, h is the time-independent solution obtained for the h field and p << h. Similarly, 
the time-independent fields / and g in equation (131!) are written as / and g. Substituting 
equation (l32l) in (|3T1) and linearizing the equation for p(r), we get 

u;^p=-[f^^ + l±]p+ [3K^~h^ + k{k + X)f - kX,P] p (33) 

ar'^ r ar 

Equation fl33|) has the form of a one-dimensional Schrodinger equation with a potential u{r) 
given by 

u{r) = 3K^h{r)^ + k{k + \yg{rf - nXifirf (34) 

Looking for imaginary modes of the frequency uo now reduces to finding negative-energy 
bound states for this potential. Figure [5] depicts the behavior of the potential for values of 
couplings which lie in both the stable and unstable regions for |V"j^^* ""-sj^ Referring to table 
[H when K = 1.3 and Ai = 0.65, \Vi ""^•') is stable for A = 1.5. This is in agreement with 
the fact that the minimum of the potential for A = 1.5 has positive energy. As A reduces 
towards A = 1.43, notice that the minimum of the energy starts reducing until it starts 
becoming negative near A = 1.43. This once again confirms the numerical result that vacua 
of type \Vi ^^"'^ ) become unavailable for values of A < 1.43. We have similarly computed 
the potential u{r) for critical values of n and A, and found similar solutions confirming the 
results stated in table [B 

Returning to our comment on the early Universe setting at the end of sec. lIVBt the field 
/ (the VEV of P and A^) has negative effective mass-squared, M2 < 0. High temperature 
corrections from gauge sector add a term of the form AT'^ with A > to this effective 



mass-squared [37|, l38|, which would drive the temperature dependent value /r to zero at 
high temperature, restoring the gauge symmetry. Below the symmetry breaking critical 
temperature /t becomes non-zero. As long as this value remains small, the effective mass- 
squared for the h in eq. fl30|) remains positive. This makes |V^i) a valid local minimum and 
admits a vortex sector \Vi "^^'y At an even lower temperature Jt grows sufficiently large, 
and the roll-over transition sets in. 

VI. CONCLUSION 

One of the main requirements in the minimal models of gauge mediated supersymme- 



try breaking is that the messenger squarks do not develop VEVs |22[]. In this paper, we 
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FIG. 5: The equivalent Schrodinger potential of equation (|33p with k = 1.3, Ai = 0.65, and different 
values of A. A positive energy minimum results in a stable solution. When a negative energy bound 
state is possible, vacua of type \Vi "^^'"^ become unavailable. 

have shown an example in which cosmic strings are present in the metastable vacua. These 
configurations can roll-over classically into undesirable configurations which asymptote into 
vacua in which (q) = (q) ^ 0. This happens for certain ranges of the couplings k, A and Ai, 
including the phenomeno logically interesting region in which Ai < 0.1. For a general mes- 
senger group Gm stable cosmic strings would exist depending on the choice of representation 
of the fields P and N. Furthermore, we have evidence that if Gm is global the instability issue 
becomes more severe. In all such models, it should be possible to obtain new constraints on 
the coupling constants from detailed simulations of the cosmic strings. Such studies provide 
tighter constraints on the parameters as compared to those from quantum tunneling effects. 
In principle, these techniques are applicable to the SUSY breaking sector itself provided 
the model supports solitonic solutions. This aspect of the problem will be dealt with in 
future works. 
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